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Abstract. We give a simple, but accurate method that 
can be used to account for illumination in compact bi- 
nary systems which have a low-mass companion, even if 
spherically symmetric illumination of the secondary star 
(not necessarily on the main sequence) is not assumed. 
This is done by introducing a multiplicative factor $ in 
the Stefan-Boltzmann surface boundary condition, which 
accounts for the blocking of the intrinsic secondary flux 
by X-ray heating of the photospheric layers. Numerical 
fits and tables for <!> are given for unperturbed effective 
temperatures in the range 2500 - 5600 K and logg in the 
range 1.0 - 5.0 
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1. Introduction 

The effects of irradiation of the low mass companion in a 
low-mass X-ray binary or a cataclysmic variable have re- 
ceived much attention in the recent years. In cataclysmic 
variables, the external illumination flux due to accretion 
onto the primary compact object is comparable to the in- 
trinsic stellar flux produced by nuclear reactions, whereas 
this external flux exceeds by orders of magnitude that of 
the secondary in low-mass X-ray binaries. A fraction of 
the illumination flux is absorbed in optically thin regions 
of the secondary photosphere, and this may result in the 
formation of a wind (see e.g. London et al. 1981; Ruder- 
man et al. 1989; Tavani & London 1993), which can in turn 
affect the long term evolution of the system. We shall not 
discuss this effect here, but instead consider energy depo- 
sition below the photosphere of the secondary. 

Early calculations in the case of low-mass X-ray bina- 
ries (Podsiadlowski 1991; Harpaz & Rappaport 1991) in 
which it was assumed that the secondary is illuminated in 
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a spherically symmetric manner showed a quite dramatic 
effect. The intrinsic luminosity of the secondary star is 
very efficiently blocked by the impinging radiation in the 
outer convective layers, so that the secondary must ex- 
pand until it becomes almost fully radiative. The secular 
evolution is then quite different from the standard case, 
namely the systems rapidly evolve towards longer periods, 
with mass transfer rates close to or above the Eddington 
limit. This would have accounted for the fact that the pe- 
riod distribution of LMXBs seems to be significantly dif- 
ferent from that of cataclysmic variables, showing a lack of 
systems at short periods. This would have also increased 
by a large factor the number of LMXBs, making it com- 
patible with what is required to account for the observed 
density of millisecond pulsars which are believed to be the 
descendants of LMXBs (Frank et al. 1992). However, the 
spherically symmetric assumption is obviously incorrect, 
and it was later shown by Hameury et al. (1993) and con- 
firmed by Harpaz & Rappaport (1995) that the unillumi- 
nated side can efficiently cool the secondary. The secular 
evolution of LMXBs is not drastically different from the 
unilluminated case, although the short term behaviour is 
very significantly affected by illumination, and exhibits on 
and off states, with relatively high mass transfer rates. 

In the simpler context of cataclysmic variables, Rit- 
ter et al. (1995, 1996a) showed that irradiation-induced 
mass transfer cycles could also be present in these sys- 
tems; King (1995) and King et al. (1995, 1996) discussed 
from a very general point of view the occurrence of such 
mass transfer cycles in CVs. The existence of these cycles 
could be responsible for the observed spreading of mass 
transfer rates M for a given orbital period P, whereas 
models would predict a fair correlation between M and 
P. Ritter et al. (1996b) showed that, in the bi-polytrope 
approximation (Kolb & Ritter 1992), the stability crite- 
rion depends critically upon the variation of the effective 
temperature of the illuminated star as a function of the 
irradiating flux. For modeling the response of the stellar 
surface to illumination, they used a very simple one zone 
model for the superadiabatic layers of the low-mass sec- 
ondary. Here, we use detailed stellar models to calculate 
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the response of the secondary, namely the secondary lu- 
minosity as a function of irradiation, for a range of values 
of the surface gravity and unperturbed effective temper- 
ature. We give our results in a tabular form that can be 
used to calculate the evolution of a compact system in the 
presence of illumination. The advantage of this formula- 
tion is that it can be used in bi-polytropic codes, which 
are much faster that full stellar codes, and thus allow the 
exploration of a wide range of parameters, but have to be 
calibrated. 



2. Reaction of the stellar surface to external 
irradiation 

The influence of heating of subphotospheric layers by an 
external irradiation flux is entirely described by a mod- 
ification of the boundary condition, as all of the ab- 
sorbed X-ray flux is thcrmalized before being re-emitted. 
In the plane parallel approximation, the standard Stefan- 
Boltzmann law L = AuR^crT^ where L is the surface 
luminosity of the secondary star in the absence of illu- 
mination (note that L need not be equal to the nuclear 
luminosity), R 2 its radius, and T c ff j0 is the effective tem- 
perature of the unilluminated star has to be replaced by 
(Ritter ct al. 1996b): 
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where x = F- m /aT^ s is the normalized illuminating flux, 
and 
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Here T c g is the effective temperature of a stellar surface 
element subject to the illuminating flux x. The function 
G(x), in the range 0-1, describes the blocking of the in- 
trinsic stellar luminosity as a result of illumination. F; rr in- 
cludes only the fraction of the flux that is deposited below 
the photosphere; as mentioned above, energy deposition in 
optically thin regions does not result in modifications of 
the internal structure of the star. 

In the case of low-mass secondaries which have a con- 
vective envelope, the entropy deep in the envelope remains 
constant over the whole surface of the star, and varies 
slowly with time, on the Kclvin-Hclmholtz time of the 
whole envelope. By contrast, the superadiabatic and radia- 
tive outer layers have a very short thermal time, and are 
in thermal equilibrium. Moreover, Gontikakis & Hameury 
(1993) have shown that the flux, in the plane parallel ap- 
proximation, docs not vary with depth in the convective 
zone. The problem of determining the reaction of the star 
to illumination thus reduces to finding the structure of an 
initially convective layer in thermal equilibrium when one 
changes the outer boundary condition while keeping the 
entropy at its base constant. This is easily done solving the 
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Fig. 1. Function G(x) for secondaries on the main sequence, 
with a mass 0.8 (a), 0.6 (b), 0.4 (c) and 0.2 M (d) 

standard stellar structure equations, which, in the plane 
parallel approximation write: 



d\ogP 
dlogS 
dlogT 
dlogP 
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(3) 



where P is the pressure, T the temperature, V the tem- 
perature gradient calculated using the mixing length ap- 
proximation in the convective zone, £ the column density, 
and g the surface gravity, assumed to be constant. The 
energy flux F is also assumed to be constant throughout 
the layer. The layer is integrated down to a depth S max of 
10 6 g cm -2 , which is deep enough that the departure from 
adiabaticity is negligible; our results have been found to be 
independent of the particular value of S max - The surface 
boundary conditions are standard: 



kP = 2/3.g, F + F- m = aT A 



(4) 



In the absence of illumination, these equations are inte- 
grated for a given g and T e g, and give the entropy So 
deep in the convective zone. In the presence of an irradi- 
ating flux, one adds the condition S = So at the base of 
this layer. The set of equations (3) are integrated using 
the method, equation of state and opacities described in 
Hameury (1991). This gives F, and thus G as a function 
of g, T e fffi, and x. Figure 1 shows G(x) for main sequence 
stars of various masses. A significant difference can be 
seen between very low mass secondaries (less than 0.4 M Q ) 
and more massive ones. This difference is due to the low 
value of the energy flux that has to be carried through the 
convective zone, and hence to a small deviation to adia- 
baticity, even in the outer superadiabatic layers. For those 
low-mass stars, G(x) is not very different from the relation 
G(x) — max(l — x, 0) that one would obtain assuming that 
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Fig. 2. Functions g(x) = —dG(x)/dx and xg(x) for secondaries 
on the main sequence with masses 0.8 (a), 0.6 (b), 0.4 (c) and 
0.2 M e (d). 

there is no superadiabatic region, so that the temperature 
gradient and thus the surface temperature is fixed. 

The stability analysis of mass transfer involves the 
derivative of the function G(x) (Ritter et al. 1996c), 
g{x) = —dG(x)/dx. More precisely, the criterion for sta- 
bility against irradiation-induced mass transfer is: 



C = 2(( s -( R )?^T- 1 >2 S xg(x) 



(5) 



M 



where ( s is the adiabatic mass radius exponent of the sec- 
ondary, tkh,o the Kelvin-Helmholtz time, the mass 
transfer time scale, T a dimcnsionclcss function which 
scales roughly as the inverse of the mass of the secondary 
convective envelope, s ~ 1/2 the fraction of the secondary 
exposed to illumination, and £r is the mass radius expo- 
nent of the Roche radius. This criterion is deduced from 
the assumption that the mass-radius exponent of the sec- 
ondary under the sole influence of illumination must be 
less than ( s — (r; Eq. (5) means that, for instability to 
occur, one must satisfy both conditions of sufficient illu- 
minating fluxes and sensitivity of the secondary response 



to irradiation, i.e. large value of g(x). C does not depend 
on the very detailed structure of the illuminated star, and 
is easily accessible in the bi-polytropic approximation. The 
response of the secondary is entirely contained in the func- 
tion xg(x), which is plotted in Fig. 2. The small discon- 
tinuities are due to linear interpolations in determining 
the opacity. A remarkable characteristic of this function 
is that (1) g(x) is a monotonously decreasing function , 
which is less than unity, and (2) that xg(x) has a max- 
imum of the order of 0.5, whatever the secondary mass, 
although the position of this maximum does depend on it. 
Thus, there will be no cycles if C > s (note that this is a 
sufficient but not necessary condition for stability). 

It must be stressed that this procedure is valid only 
for stars which have a convective envelope, even though 
the whole star can be far out of thermal equilibrium. This 
would in fact be the main limitation in following the evo- 
lution of systems in which the effect of irradiation is so 
strong that at some point the star becomes fully radiative. 
This was a natural outcome of models assuming spheri- 
cally symmetric illumination, but Hameury et al. (1991) 
have shown that this is no longer the case when one ac- 
counts for the asymmetry of irradiation. 

These calculations can easily be generalized to stars 
which arc not on the main sequence, provided they still 
have a convective envelope. Given the effective tempera- 
ture T e ff and surface gravity g, the set of equations (3) can 
be integrated in the absence of illumination and provide 
the unperturbed outer structure; the same procedure as 
previously is then applied. Table 1 (available only in elec- 
tronic form at the CDS via anonymous ftp 130.79.128.5) 
gives — log 10 (G) (column 4) as a function of log g (column 
1), logTcff (column 2) and log(l + x) (column 3). 

If one assumes that the irradiation source is a point lo- 
cated at the distance a from the center of the secondary, 
Eq. (1) can be easily integrated over the whole stellar sur- 
face, which gives, following the notations of Ritter ct al. 
(1996b): 



L = 2^i? 2 2 aT c % 
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G(x(9))sin9d9 



,(6) 



where 9 is the colatitude with respect to the substellar 
point of a point on the secondary surface, /2(g) = R2/CI, 
with q = M2/M1 the mass ratio, and the normalized irra- 
diation flux x is now given by: 



x(9) = x h(9) 
with 

r] GM\M 
4^ R x a 2 



x 



- 



(7) 



(8) 



where r] is an efficiency factor accounting for the albedo of 
the secondary, the fraction of energy deposited in optically 
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Fig. 3. Evolution of a cataclysmic variable using either a two 
hemispheres model with standard boundary conditions (top 
panel) , or a spherically symmetric model with condition (9) 

thin regions or the anisotropy of emission at the surface 
of the primary, and 



h{6) 



cos 6 — fi 



(l-2/ 2 cos0 + /f) 3 / 2 • 
Equation (6) can thus be written in the form 
L = 47ri?2aT c 4 ff $(x ) . 



(9) 



(10) 



We have integrated numerically Eq. (6), and the results 
can be fitted within 1 % by: 



y = ln(l + x ) 

n^o) = — ^— + -^r~ 



exp(aij/ + a 2 y 2 + a 3 y 3 ) 



Table 2 (available only in electronic form at the CDS via 
anonymous ftp 130.79.128.5) gives the coefficients m, a 2 , 
and 03 as a function of logT o ff, and \ogg. For arbitrary 



values of logr e ff j0 and log g, a linear interpolation of y 
(not of the coefficients a±, a 2 and 03) can be done. 

The stability criterion now involves the integral of the 
function xg(x); mass transfer will be stable if: 



C > max 
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We have calculated this maximum using the values of G(x) 
as determined above, and found that this maximum is 
again almost independent of the secondary mass, and is 
smaller by about a factor 2 than the value 2s found assum- 
ing uniform illumination over a fraction s of the secondary, 
being equal to (0.2 - 0.3) (1 - / 2 ). 

We have compared evolutionary calculations using ei- 
ther a spherically symmetric code with a boundary con- 
dition given by (10), or the code described in Hameury 
et al. (1993) in which both hemispheres with a coupling 
term are modeled. The results are given in Fig. 3, for a 
cataclysmic variable with Mi = 1 M and M 2 — 0.8 M Q ; 
the white dwarf radius is 5 x 10 8 cm, and the efficiency r] 
is 0.1 (spherically symmetric model) and 0.09 (two hemi- 
spheres model). Both are very similar, i.e. the mass trans- 
fer rate exhibits damped oscillations; the different values 
of r\ were chosen so as to obtain the same initial maximum 
value of the mass transfer rate. The differences (slightly 
shorter damping time and 77) are mainly due to the fact 
that in the two hemispheres model, the illuminating flux is 
assumed to be constant over the heated regions, whereas 
the boundary condition (10) includes an angular variation 
of this flux, which has a strongly non-linear effect. 

For X-ray luminosities close to the Eddington limit, 
wc predict that the illumination effect is small, as xg(x) 
is small for large values of x, and as a significant fraction 
of the secondary surface might be shielded from illumina- 
tion by the accretion disc (Ritter et al., 1996b). However, 
Eq. (10) might significantly underestimate the effect of 
heating, as the heat flux in the illuminated regions can be 
negative, and be of the same order of magnitude as the in- 
trinsic stellar flux. The unilluminated layers have then to 
re-radiate this additional flux, which may, in some cases, 
be sufficient to also block the intrinsic stellar flux. This 
is equivalent to having x ~ 1 in these regions which then 
become very sensitive to illumination. This effect, respon- 
sible for the short outbursts of mass transfer obtained by 
Hameury et al. (1993), is indeed not accounted for by the 
boundary condition (10), but depends sensitively on the 
strength of the coupling between both hemisphere, i.e. on 
the circulation timescale which is poorly known. 



(11) 3. Conclusion 



The effects of X-ray illumination of the secondary star in 
a CV or a LMXB can be accounted for by a modification 
of the Stefan-Boltzmann law which we give in a tabular 
form. This approximation is quite accurate in CVs, but 
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may lead to a significant underestimate of the effects of 
illumination in LMXBs, when the luminosity is close to 
the Eddington limit. This modified boundary condition 
can be used in bi-polytropic codes, which are much faster 
than full stellar codes, and hence allow the exploration of 
a larger parameter space. 
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